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Maximally decimated filter banks have been extensively studied in the past. A 
filter bank is said to be under-decimated if the number of channels is more than the 
decimation ratio in the subbands. A maximally decimated filter bank is well known 
for its application in subband coding. Another application of maximally decimated 
filter banks is in block filtering. Convolution through block filtering has the advan- 
tages that parallelism is increased and data are processed at a lower rate. However, 
the computational complexity is comparable to that of direct convolution. More 
recently, another type of filter bank convolver has been developed. In this scheme, 
the convolution is performed in the subbands. Quantization and bit allocation of 
subband signals are based on signal variance, as in subband coding. Consequently, 
for a fixed rate, the result of convolution is more accurate than is direct convolu- 
tion. This type of filter bank convolver also enjoys the advantages of block filtering, 
parallelism, and a lower working rate. Nevertheless, like block filtering, there is no 
computational saving. 

In this article, under- decimated systems are introduced to solve the problem. 
The new system is decimated only by half the number of channels. Two types 
of filter banks can be used in the under- decimated system: the discrete Fourier 
transform (DFT) filter banks and the cosine modulated filter banks. They are well 
known for their low complexity. In both cases, the system is approximately alias 
free, and the overall response is equivalent to a tunable multilevel filter. Properties 
of the DFT filter banks and the cosine modulated filter banks can be exploited 
to simultaneously achieve parallelism, computational saving, and a lower working 
rate. Furthermore, for both systems, the implementation cost of the analysis or syn- 
thesis bank is comparable to that of one prototype filter plus some low-complexity 
modulation matrices. The individual analysis and synthesis Biters have complex co- 
efficients in the DFT Biter banks but have real coefBcients in the cosine modulated 
Biter banks. 


I. Introduction 

The M channel maximally decimated filter bank shown in Fig. 1 has been studied extensively in 
[1-8]. A filter bank is said to be under-decimated if the number of channels is more than the decimation 
ratio in the subbands. When the system in Fig. 1 is alias free, it is a linear time invariant system with 
transfer function T(z), as indicated in Fig. 1. In the following discussion, T(z) will be called the distortion 
function, or the overall response. 

A maximally decimated filter bank is well known for its application in subband coding. Another 
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of the synthesis bank need not be the same. To be more specific, let P 0 (z) be the prototype filter of 
the analysis bank and Q 0 (z ) be the prototype filter of the synthesis bank. The filters P k (z) and Q k (z), 
k = 1, 2, • • • , 2M - 1, are, respectively, the shifted versions of P 0 (z) and Q 0 (z). 


P k (z) = P 0 (zW k ), Q k (z ) = Q 0 (zW k ), k = -M, -M + 1, • ■ . , M - 1 


Notice that unit circle P k (z) is only a shift of P 0 (z) by kn/M, since P k (e juj ) = P 0 (e :? ' (a,_fc7r/M) ). Figure 3 
shows this relationship. The analysis filters and synthesis filters of the DFT filter bank have the following 
form: 


H k {z) = a k P 0 (zW k ), and F k (z) = a* k Q 0 (zW k ), W = 


IP 0 (e^)f I Pi(e/")l 


\P k (efc)\ 


-n/2M 0 n/2M n/M kn/M 

Fig. 3. Magnitude responses of P k {z). 




(1) 


(The definition of a DFT filter bank here is slightly different from the conventional DFT filter banks [7].) 
It follows that H k (e JU> ) is just a shift of P 0 (e^) by kn/M except for a scalar. The same holds for the 
synthesis filters. 


We now show that with proper design of the two prototypes, this DFT filter is approximately alias 
free and the overall response is equivalent to a tunable multilevel filter. Moreover, the overall response 
can be a real-coefficient linear-phase filter as desired. Efficient implementation of the DFT filter bank 
will also be discussed. 

A. Suppression of Aliasing Error 

Consider the under-decimated system in Fig. 2, a 2M-channel filter bank with decimation ratio M. 
The suppression of aliasing error due to downsampling in the subbands can be explained pictorially. Take 
the first subband as an example. Because of decimation followed by expansion, there will be M — 1 image 
copies of H 0 (z), as shown in Fig. 4. We can see from Fig. 4 that these image copies will be suppressed 
if both H 0 (z) and F 0 (z ) have stop-band edges less than n/M. When the spectral supports of F 0 (z) and 
the image copy of H 0 (z ) do not overlap, the aliasing error will be suppressed to the level of the stop-band 
attenuation of H 0 (z ) or F 0 (z), which is equivalent to the stop-band attenuation of P 0 (z ) or Q 0 (z). In the 
other subbands, the same reasoning for aliasing suppression applies. 

We now present the mathematical counterpart of the above discussion. The output X(z) is related to 
the input X (z) by 


M — l 

X(z) = Ai(z)X(zW 2i ) 

i~0 


( 2 ) 


The alias transfer function, Ai(z), is defined as 
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ltf 0 {e/®W- 4 )| \H 0 {e^W- 2 )\ 


IH 0 (e^°W 2 )l 

IH 0 (e/ co W 4 )l 


0 

2k /M 



IF 0 (e/“)l 



0 


Fig. 4. Image copies of Hq ( z ) due to decimation followed by expansion and the spectral support of Fq ( z ). 

^ 2M — 1 

M*) = H k (zW 2i )F k (z) ( 3 ) 

k = 0 


The system in Fig. 2 is alias free if Ai(z) — 0 for i — 1 , 2, * • • , M — 1 . With analysis filters and synthesis 
filters chosen as in Eq. (1), Ai(z) can be written as 

1 2M— 1 

Ai ^ = M £ M 2 Po(zW 2i+k )Q a (zW k ) ( 4 ) 


Assume the nonadjacent bands of Pq(z) and Q 0 (z ) do not overlap, i.e., 

-Po(e 3 "VK 2i )Oo(0«0, i = (5) 

This assumption is reasonable if Po(z) and Qo(z) have stop-band edges less than 7 r/M and large enough 
stop-band attenuation. Equation (5) gives us 

fc = 0, 1, • ■ • , 2M — 1, z = 1, 2, • • • , M — 1 

which implies Ai(z ) w 0,z ^ 0. We conclude that the DFT filter bank is almost alias free. Also notice 
that the degree of alias suppression improves with the stop-band attenuation of the two prototypes. 

B. The Overall Response of the DFT Filter Bank 

For a 2 M- channel system decimated by M as shown in Fig. 2, the distortion function T(z) can be 
expressed as [7] 

^ 2M— 1 

T (*) = ^ £ H k (z)Fk(z) ( 6 ) 

A:=0 


Let Rq(z) = P 0 (z)Q 0 (z). Substitute the expression of H k (z) and F k (z) in Eq. (1); then 

1 2M-1 

= M £ ( 7 ) 

k = 0 
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When i2o(z) is a Nyquist (2 M) filter, it can be shown the addition of |a fc | 2 tf 0 (e i(w-A:7r/M )) in Eq. (7) will 
not result in any bumps or dips in the response of T(z ) because of the Nyquist property of Rq{z). The 
definition of a Nyquist filter is given in Section I. A detailed explanation can be found in [7]. 


With Eq. (7), we can plot a typical magnitude response of T(z), as in Fig. 5, which shows that the 
overall response is equivalent to a multilevel filter. Since the value of a k can be chosen freely, T(z) 
is actually a tunable multilevel filter. In particular, we can get lowpass filters with stop-band edges 
adjustable in integer multiples of 7 r/M. 



Fig. 5. A typical magnitude response of T(z), a multilevel filter. 


Remarks: 


(1) If Rq(z) is a real-coefficient filter and we choose a k = o, 2 M—k to he some real number, k = 
1, 2, • • • , M, it can be verified that the resulting T(z) is also real- coefficient. 

(2) Let Ro(z) be linear phase with order N r , a multiple of M. In this case, R 0 (zW k ) is linear phase. 
By Eq. (7), this implies that T(z ) has linear phase. 

Summarizing, we have shown that if R 0 (z ) is Nyquist (2 M) and Eq. (5) is valid, the DFT filter in 
Fig. 2 is nearly alias free and the overall response is equivalent to a tunable multilevel filter. 

C. Implementation of the DFT Filter Bank 

There exists an efficient implementation for the DFT filter bank. To see this, express the prototype 
Pq(z) in polyphase representation, i.e., 


2M— 1 

Po(z) = £ E i (z 2M )z- i ( 8 ) 

i=0 

where Ei(z) is the zth type 1 polyphase component of P 0 (z ) [7]. The analysis filters can be rewritten as 

2M-1 

Hk(z) = a k P 0 (zW k ) = a k ^ Ei(z 2M )W- ki z-\ k = 0, 1, • • • , 2M - 1 (9) 

7=0 

Let 


h(z) — [Hq(z) Hi(z) • • • H2m-i{z)] t 


( 10 ) 
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The matrix representation of Eq. (9) is 


(a 0 

0 • 

0 

\ 

/ Eo(z 2M ) 0 

0 x 

0 

h (*) = 


ai • 

0 


w* 

0 Ei(z 2M ) • 


V 0 

0 * 

• • «2M- 

J 

1 

\ 0 0 

•• E2M-l(z 2M ) / 


&2m(z) ( 11 ) 


Observing Eq. (11), we can draw the polyphase implementation of the analysis bank as in Fig. 6. The 
implementation cost is that of the prototype filter P 0 (z) plus a DFT matrix. The same holds for the 
synthesis bank. The computational complexity of the analysis bank is comparable to that of the analysis 
prototype filter plus one DFT matrix. Notice that all the computations involved in the filter bank are 
performed after the M-fold decimators; lower rate and lower complexity are achieved at the same time. 



Fig. 6. Efficient implementation of the analysis bank of the 2/lf-channel DFT filter bank. The DFT matrix, W, 

is of size 2M x 2M. 


III. Cosine Modulated Filter Banks and Application in Tunable Multilevel Filtering 

In the DFT filter bank described in the previous section, the analysis and synthesis filters have complex 
coefficients. If it is desirable for the individual filters to have real coefficients, then we can use the new 
under-decimated cosine modulated filter bank to be discussed in this section. 

The system in Fig. 2 is said to be a cosine modulated filter bank if all analysis and synthesis filters 
are generated by cosine or sine modulation of one or two prototype filters. In this section, we introduce 
two classes of new under-decimated cosine modulated filter banks. The systems are nearly alias free. 
Aliasing error decreases as the stop-band attenuation of the prototype increases. Individual analysis and 
synthesis filters have real coefficients. We can design the prototypes so that the overall response of the 
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filter bank is a linear-phase multilevel filter. Furthermore, there exists efficient implementation of this 
cosine modulated filter bank. The implementation cost of the analysis bank is that of the prototype filter 
plus two DCT matrices (Appendix). The complexity of an M x M DCT matrix is only of the order 
Aflog(M) [11]. The same holds true for the synthesis bank. Two types of filter stacking can be applied 
in the new under-decimated cosine modulated filter bank. In our discussion, a cosine modulated filter 
bank will be identified as type 1 or type 2 according to the stacking of its analysis and synthesis filters. 

A. Type 1 Cosine Modulated Filter Bank 

1. Construction of the Type 1 Cosine Modulated Filter Bank. In the cosine modulated filter 
bank, all analysis and synthesis filters have real coefficients. Each filter has positive and negative spectral 
occupancy as opposed to single-sided spectral occupancy in the DFT filter bank. This incurs a problem 
that we do not have in the DFT filter bank. Details and a proposed solution of this new problem will be 
given below. 

Let P Q (z) and Q 0 (z ) be respectively the prototype filters of the analysis bank and the synthesis bank. 
The definitions of P k (z ) and Q k (z ) are as in Section II. To get real-coefficient analysis and synthesis filters 
from the prototypes, we can combine Pk(z ) and P- k (z): 


H k {z) = a k P k (z) + a* k P_ k (z ), F k {z) = b k Q k (z ) + b* k Q_ k (z), k = 1 , 2 , • ■ • , M — 1 

Since Pq(z) and Pm{z ) are already real filters, we can directly choose 


H k {z) = 2 a k P k (z), F k (z ) = 2 b k Q k (z), k = 0 or M 


Figure 7 shows the spectral supports of the analysis filters for the type 1 cosine modulated filter bank. 
The stacking of the spectral supports of the synthesis filters is similar. 


h M h M- 1 H i H o H 

\ 

• • • 


%-1 h M 


-Jt -(M-1)jt/yW -n/M 0 n/M (M-1 )k/M n 

Fig. 7. Normalized magnitude responses of the analysis filters of the type 1 cosine modulated filter bank. 


Aliasing error due to M-fold decimation in the Oth and the Afth subband can be suppressed on the 
synthesis side as we did in the DFT filter bank. The situation in the other subbands is different because 
now H k {z) and F k (z ), k — 1, 2, • • ■ , M — 1, are bandpass filters. Referring to Fig. 8, decimation by M 
followed by expansion by M in the subbands will cause one image copy of P k (z) to overlap completely 
with Q_ k (z), k = 1 , 2, • • • , M — 1 . This type of aliasing error cannot be suppressed in the synthesis bank. 

Our solution to this problem is to introduce a second subsystem that has exactly the same aliasing 
error to cancel the existing one. Let the second subsystem have analysis filters H {z) and synthesis filters 
F k ( z )i k = 1> 2, • • • , M — 1. To create the same aliasing error, the filters of the second subsystem are 
required to have similar stacking as that of the first subsystem. In particular, 


H 'k( z ) = a 'k p k(z) + 4*P_ fc (z), F' k {z) - b' k Q k {z) + b k Q_ k (z), k = 1, 2, • • • , M — 1 


The configuration of the analysis filters in the second subsystem is shown in Fig. 9. Notice that the 
spectral occupancy of Hq^z) and Hm{z) are not covered in the second subsystem. 
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(a) 

I P k (e y ® W" 2/f - 2 )l I P k (e/® W~ 2k )\ I P k (e^® V\r 2 * +2 )| I P k (e'® V\r 2 )l \P k ( e l® W 2 )l 



Fig, 8. Type 1 cosine modulated filter bank: (a) image copies of P^z) due to decimation followed by 
expansion and (b) the spectral support of F^z). 



Fig. 9. Normalized magnitude responses of the analysis filters of the second subsystem 
in the type 1 cosine modulated filter bank. 


The setup of the new system is now complete and is shown in Fig. 10. It is a connection of two 
subsystems. The first subsystem has M + 1 channels, and the second subsystem has M - 1 channels. 
The whole system is under-decimated; it has 2 M channels but is decimated only by M. The analysis and 
synthesis filters can be summarized as follows: 



Fig. 10. The setup for the type 1 under-decimated cosine modulated filter bank. 
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H k (z) = 2 a k P k (z) 


k = 0,M 


H k (z)=a k P k (z) + a* k P- k (z), 

k- 1,2,-. 

• , M — 1 

H' k (z)=a' k P k (z) + a' k *P_ k (z), 

k = 1,2,*- 

• , M — 1 

F k (z) = 2 b k Q k {z), 

k = 0,M 


Fk(z) — b k Q k (z) + b k Q- k (z), 

k = 1,2,-- 

■ ,M — 1 

Fk( z ) — b' k Q k (z ) + b k Q- k (z ), 

k = 1,2,-- 

• , M — 1 


The values of a k , a' k , b k , and b' k will be determined later. 

In the following we show that with proper design of the prototypes and appropriate choices of a k , a' k , 
b k , and b k , this filter bank is almost alias free. The overall response of the type 1 cosine modulated filter 
bank can be designed to be a linear-phase tunable multilevel filter. 

2. Cancellation and Suppression of Aliasing Error. As we mentioned in the construction of 
filters, the aliasing error in the Oth and the Mth subbands will be suppressed in the synthesis bank. Only 
the subbands with bandpass filters require alias cancellation. The physical picture is as follows: 

Consider the fcth subband, 1 < k < ‘M - 1. Due to decimation followed by expansion, P k {z) has M - 1 
image copies and P^ k (^z) also has M — 1 image copies. The image copies of P k [z) will be suppressed 
by Qk(z ) provided that both Pq(z ) and Qo(z ) have stop-band edges less than tt/M and large enough 
stop-band attenuation. Of the M — 1 image copies of P k (z), M — 2 of them are in the stop band of 
Q~ k (z) and, hence, will be suppressed by Q- k (z), as depicted in Fig. 8. However, one of the image copies 
of P k (z ) will fall on top of the spectral support of Q~ k (z). Only this image copy requires cancellation. 
In the fcth subband of the second subsystem, the same aliasing occurs. It can be shown that the aliasing 
error of the second subsystem cancels that of the first subsystem if the values of a k , b k , a' k , and b' k are 
chosen properly. Mathematical proof of this claim is as follows. 

With filters constructed as in Eq. (12) and the expression of alias transfer functions in Eq. (3), we 
have 

= jj K’M + 4 2) M + A f\z) + A?\z)) (13) 

where 

M— 1 

A^iz) = (a 0 b 0 + a* 0 b 0 ) P 0 (zW 2i ) Q q (z) + (akb k P a' k b' k ) P k (zW 2i ) Q k (z) 

k = 1 

M - 1 

z) = J2 (“A + *ZK) P-k (^) Q k (z) 

k= 1 

M — 1 

= E + a X) Pk (^ 2i ) Q-k(z) 

k= 1 
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4 4) W = K»5 + “oK) Po (zW 2i ) Qo(z) + J2 M K + 0,'X) P-k (zW 2i ) Q-k(z) 


M - 1 


A:=l 


Assume P 0 (z ) and Q 0 (z) satisfy Eq. (5). It follows that A { p{z) « 0 and Af\z) « 0. Let a k , k = 0, * • ■ , M 
be real and choose 


bk — ak, k = 0, • ■ • , M 

a k = -J a k, b' k = ja k , * = 1,---,M-1 


(14) 


If a*, 6*,a£, and 6J. are chosen as above, it can be verified that a* k b k + of 6J. = 0, which implies Af } fy) = 

(z) = 0. So the condition A f (z) « 0 is ensured, provided that Eq. (5) is valid and a k ,b k ,a' k , and bi 
are chosen according to Eq. (14). 

With Eq. (14), we can write down the time domain description of the analysis and synthesis filters for 
the type 1 cosine modulated filter bank. Let p Q (n) be the impulse response of P 0 (z) and q 0 (n) be the 
impulse response of Qo(z). 


h k (n) — 2a k p 0 (n)cos ( knir/M ), 
h' k (z) = 2akPo(n)sm ( knix/M ), 
fk(z) = 2akqo(n)cos ( krm/M ), 
f k (n) = 2a,kqo(n) sin (krnr/M), 


fc = 0,l,-,M 

k = 1,2, • ■ • , M — 1 

* = 0,1, ---.Af 

* = 1, 2, ■ ■ ■ ,M — 1 


(15) 


From the expression in Eq. (15), we can see that each individual filter is a sine or cosine modulation of 
the prototype filters. 

3. The Magnitude Response of the Overall Response T{z). Using Eqs. (12) and (6), we get 


T (e JW ) 


2 

M 


M 


M 2 ((a 0 (e^ ( " +far/M )) + Bo 

fc— 0 ' ' ' ' 


(16) 


The above expression for the overall response is similar to that in the case of DFT filter bank, Eq. (7). If 
Ro{z) is a Nyquist (2 M) filter, this is a tunable multilevel filter bank like in a DFT filter bank. 

4. The Phase of the Overall Response T(z). The overall response T(z) has linear phase provided 
that Rq{z) is linear phase and N r , the order of Rq(z), is an even multiple of M. The reason is given 
below. The linear phase property of Rq(z) entitles us to write 


R 0 (e JW ) = e~^ N ^ 2 R{u) 

where R(u) is a real- valued function. Substitute Eq. (17) into Eq. (16) and we get 


(17) 
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T (e jtv ) 


A M / / 

±e ~^' 2 Kl 2 ( ( w ~ k1 ' 

t=o ' \ M 


q_ e ~jknN r /2M ft 


f CJ + kir \ \ 

V M )) 


which has linear phase when N r is a multiple of 2 M. 

f Notice that if Q 0 (z) is the time reversed version of P 0 (z), i.e., Q 0 (z) = z~ N r>Po(z), then F k (z) and 
F k(z ) are the time-reversed version of H k (z) and H' k (z), respectively. With this choice of Q 0 (z), the 
overall response is 


T{z) 


Z -N p / M ____ M - 1 

-JT + E Hi(z)H' k (z) 

\k= 0 fc=l 


In this case, the overall response has linear phase regardless of the order of R 0 (z). 

5. The Phase of the Individual Analysis and Synthesis Filters. Let the prototype filter P 0 (z) 
be linear phase. If the center of symmetry of p 0 (n) coincides with that of the sine or cosine functions that 
modulate p 0 (n), then the resulting analysis filters also have linear phase. In the case of the type 1 cosine 
modulated filter bank, the condition can be further reduced. It can be verified that if N p , the order of 
p o(z), is a multiple of M, every analysis filter has linear phase. The same holds for the synthesis bank. 

When N p is an even multiple of M, all the analysis filters in the second subsystem, H' k {z),k = 
1 ) 2, • • • , M — 1 , are found to have a null at zero frequency and 7r. We can also verify that as N p is an odd 
multiple of M, the analysis filters H k (e ju ),k = 1, 2, • • • , M - 1, have zeros at zero frequency and 7 r. We 
conclude that some M — 1 of the 2 M analysis filters have zeros at zero frequency and 7 r if the analysis 
filters have linear phase. More on this property will be addressed when we present a similar result for the 
type 2 cosine modulated filter bank. 

Summarizing, we have shown that the filter bank in Fig. 10 is equivalent to a linear-phase tunable 
multilevel filter if the following two conditions hold: (1) The nonadjacent bands of P Q (z) and Q 0 (z) do not 
overlap, [Eq. (5)j, and (2) Rq(z) is linear phase and close to a Nyq.uist (2 M) filter. The implementation 
cost of the analysis bank, as will be shown in the Appendix, is the prototype filter P 0 (z) plus two DCT 
matrices. Complexity follows the buildup of the architecture; the computational cost of the analysis bank 
is that of P 0 (z) plus two DCT matrices working at an M-fold decimated rate. 

B. Type 2 Cosine Modulated Filter Bank 

In Figs. 7 and 9, we show the configuration of the analysis filter for the type 1 under-decimated cosine 
modulated filter bank. In the type 2 cosine modulated under-decimated filter bank, a different stacking of 
filters is applied. We show the new setup in Fig. 11. The filter bank can still be conceived as a connection 
of two subsystems, both with M channels. The spectral supports of the second set of analysis filters are 
exactly the same as the spectral supports of the first set of analysis filters (Fig. 12). The same holds for 
the synthesis bank. To be more specific, let P 0 (z ) and Q 0 (z ) be the two prototypes as before. Define 


Pk(z) = P Q (zW k+0 - 5 ) and Q k (z ) = Q 0 (zW k+0 b ) , W = 
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Fig. 11. The setup for the type 2 under-decimated cosine modulated filter bank. 



-n -(/W-1 )n/M -Ji /M 0 n/M ( M-\)n/M n 


Fig. 12. Normalzed magnitude responses of the analysis filters of the type 2 
cosine modulated filter bank. 

Spectral support of P k (z) is shown in Fig. 13. It is similar for Q k (z). Notice the difference between 
the definition of P k {z) and the definition of P k (z) in the type 1 cosine modulated filter bank, Fig. 3; if the 
type 1 system and the type 2 system have the same analysis prototype, then there is on the unit circle 
P k (z) a shift of P k (z ) by 7r/2M, and this comes from the extra W 0 5 in the preceding equation. 


\P 0 (e^)\ IP 1 (e/®)l IP^e/“)l 



0 n/M 2 n/M kn/M (k+^n/M 


Fig. 13. Magnitude responses of P*(z). 

We choose the analysis and synthesis filters, similar to the construction of the type 1 system, as follows: 
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H k (z)=a k P k (z) + a* k P- k (z), 

k = 0, 1, • • • ,M — 1 

Hk( z ) ~ ~3 (^ a kPk{z) ~ a kP-k( z )j ) — 0, 1, • • ■ , M — 1 


Fk{z) = 6fc§fc(^) + 6fc§_fc(^), fc = 0, 1, • • • , M - 1 

K(z) = 3 (hQk(z) - biQ- k (z )) , * = 0,l,-.-,JWf-l 


( 18 ) 


In this case, we found that the overall response is still equivalent to a tunable multilevel filter, and the 
argument of alias cancellation continues to hold after minor adjustments. More details are given below. 
The implementation cost and the computational complexity are the same as for the type 1 system. This 
can be proved in a manner similar to that used for the type 1 system in the Appendix. 

1. Cancellation and Suppression of Aliasing Error. The cancellation and suppression of alias- 
ing error is very similar to that of the type 1 cosine modulated filter bank. The image copies of Pk(z), 
Fig. 14(a), can suppressed by Qk(z), provided that P 0 (z) and Q 0 (z) have large enough stop-band atten- 
uation and their nonadjacent bands do not overlap. We can see from Fig. 14(a) that two image copies 
of Pk{z ) overlap with Q-k(z). Figure 14(b) shows the resulting aliasing that cannot be suppressed by 
the synthesis filters. As in the type 1 case, the second subsystem contributes another alias term, which 
can be used to cancel the alias from the first subsystem. This can be verified by directly substituting the 
analysis and synthesis filter in Eq. (18) into the definition of alias transfer functions in Eq. (3). 


\P k {ej (0 W- 2l< - 2 )\ \Pf^e^W~ 2k )\ \P^W- 2 )\ \P k {ej (0 W 2 )\ 



-{k+1)n/M -kn/M 0 (k-2)n/M (k+2)n/M 


FROM Q_ k 



\F k (ej<»)\ 


rV 


FROM Q k 


~(k+'\ )n/M -kn/M 


kn/M (k+1)n/M 


(b) 


~(k+-\)n/M -kn/M 


0 


co 


Fig. 14. Type 2 cosine modulated filter bank: (a) image copies of P k (z) due to decimation followed by expansion 
and the spectral support of F k (z) and (b) residual alias component that cannot be suppressed by Q k (z). 


2. The Magnitude Response of the Overall Response T{z). The overall response T(z) in the 
type 2 cosine modulated filter bank can still be expressed in a form similar to that in the type 1 case. 
Let R 0 (^) be defined as before, that is, R 0 (z) = Po(z)Q 0 (z). With analysis and synthesis filters chosen 
as in Eq. (18), the overall response is 
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T W=M E M 2R ° 

k=~M 


) 


As in the DFT filter bank, the magnitude response of T(z) will resemble that of a multilevel filter if R 0 (z ) 
is a Nyquist (2 M) filter or close to a Nyquist (2 M) filter. Again, we can tune the values of a k as desired; 
the overall response is actually a tunable multilevel filter. 

3. The Phase Response of the Overall Response T(z). In the type 1 cosine modulated filter 
bank, the phase of Tfy) is linear if Rq(z) is linear phase and N r , the order of Rq^z), is an even multiple 
of M. It is still true in the type 2 system. We can verify this by following the same steps as in the type 
1 case. 

4. The Phase Responses of the Individual Filters. Each analysis filter in the type 1 cosine 
modulated filter bank has linear phase if P 0 (z ) has linear phase and N p , the order of P 0 (z), is an even 
multiple of M. The same holds for the synthesis bank. This is still true in the type 2 case. However, 
the type 2 system is found to have the following two additional properties when the individual filters are 
linear-phase: 


(i) The filters H' Q {z) and F^z) have zeros at zero frequency, i.e., H' Q (e J '°) = F(> (e^°) = 0. 
The highpass filters H' M _ 1 (z ) and F'm- d*) have zeros at ix. 

(ii) Either H 0 (z ) or F 0 (z) has a bump at zero frequency and either H M -i(z ) or F M -i{z ) 
has a bump at n. 

Property (i): Let P 0 (z ) be linear phase and a 0 be real. Referring to the construction of filters in 
Eq. (18), the impulse responses of H k (z) and H' k {z) are, respectively, h k (n) and h' k {n ): 

h k (n) - 2a 0 p Q (n ) cos (w(k + 0.5)-^^ 

Ki n ) = 2a 0 po(n) sin (ix(k + 


Their time reversed versions are 


hk(N p - n) = 2a 0 p 0 (n) cos ( n(k + 0.5)^ - n (k + 0.5) — 

V M M 


K( N P - n ) = - 2a 0 p 0 (n)sin ^7r(A: + 0.5)-^ ~n(k + 0.5)~J 

For the above two filters to have linear phase, N r is necessarily a multiple of 2 M. If N r is an even multiple 
of 2 M, then h k {n) = h k (N p — n) and h' k {n ) = —h k (N p — n). When N r is an odd multiple of 2 M, H k (z ) 
becomes odd symmetric and H' k {z) becomes even symmetric. Without loss generality, we can assume 
that N r is an even multiple of 2 M. 

Notice that an even N p , h' k (n) = -h' k (N p - n) implies that = H' k (e *) = 0 [7], In particular, 

the lowpass filter H' 0 (z) has a notch at zero frequency, and the highpass filter has a notch at tt. 
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The other filters, H' k (z),k = 1, 2, • ■ • , M — 2, are all bandpass filters; zeros at zero frequency do not have 
a significant effect on their shapes. The same holds true for the synthesis bank. 

In the type 1 cosine modulated filter bank, we obtained similar results. In that case, when the analysis 
filters have linear phase, some M — 1 of them have zeros at zero frequency and n. But those are all 
bandpass filters; zeros at zero frequency and 7r are not of particular importance. 

Property (ii): For the flatness of the passband of T(z), Rq(z) is required to be a Nyquist (2 M) filter 
or close to a Nyquist (2 M) filter, as we have discussed previously. The Nyquist property of R 0 (z ) ensures 
that 


M-l 

e ju W k+ 0.5j = c 

k=—M 

for some constant c. Without loss of generality, we can consider c = 1. At zero frequency, only the term 
k = — 1 and the term k — 0 contribute significantly, so 

Ro (e J,r / 2M ) + flo « 1 (19) 



Since Rq(z ) is a real filter, we have 


|Ro (e^/ 2M ) | = |flo (e- 2,r/2M ) | (20) 

Combining Eqs. (19) and (20) leads to |Po (e J,r / 2M ) | > 0.5, i.e., |P 0 (e J7r / 2M ) Q 0 (e J7r / 2M ) | > 0.5. This in 
turn implies |P 0 \e jn/2M ) \ > \/oT5 or \Qq ( e ^' 7r / 2M ) | > Suppose |P 0 ( e ^ ,7r / 2M ) | > V0~5- Moreover, 

by the linear phase property of Po(^), we can write P 0 (e jui ) — e -^^ P / 2 p( w ) ; where P(w) is a real-valued 
function. With this expression, Hq (e j0 ) becomes 


Hq (e^ 0 ) = e~^ Np ^ e J' iV P 7r / 4M p £ _|_ e -jN v ir/4M p ^ ^ 

From the discussion of property (i), we know N p is an even multiple of 2 M. This gives us 

^o(e 2 °)|= 2 |P^|>^ 

This means that H 0 (z) has a bump of about 3 dB at zero frequency. If we assume \Qo (e J7r / 2M ) j > 
then F 0 (z) will have a bump of about 3 dB at zero frequency. The filters Hm-i(z) and Fm-i(z) are shifts 
of Ho(z ) and Fq(z) by 7r; if either Hq{z) or Fq(z) has a bump at zero frequency, then either Hm-i(z) or 
Fm-i{z) has a bump at zero frequency. 

In the above derivation, we assume ak is a real number. The readers will find that for complex a*,, 
properties (i) and (ii) continue to hold after minor adjustments of the above argument. 

It is noteworthy that the null of H' 0 (z ) and Fq(z) at zero frequency and the bump of H 0 (z) or Po(-z) 
at zero frequency do not affect the overall response. The overall response is a summation of responses of 
the first subsystem and the second subsystem. The bump from the first subsystem compensates for the 
null of the second subsystem to ensure the flatness of the overall response T(z). 
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IV. Design Procedures 

We now discuss four techniques for designing P Q (z) and Q 0 (z). The resulting Rq(z) = P 0 (z)Q 0 (z) is 
close to a Nyquist (2 M) filter and, at the same time, the nonadjacent bands of Po(z) and Qo(z) do not 
overlap, i.e., Pq{z ) and Qo(z) satisfy Eq. (5). The prototypes Pq(z ) and Qo(z ) designed by these methods 
can be applied to both the DFT filter banks and the two classes of cosine modulated filter banks. 

A. Eigenfilter Design 

The basic idea of this method is the following. Let P$ (z) be a lowpass filter with a stop-band edge 
less than i r/Af. (Any lowpass filter design technique can be used for designing P 0 (z).) First we find 
the condition for Qo(z ) such that the product Rq{z) = Pq(z)Qo(z) is Nyquist (2 M). Then we use an 
eigenfilter design and incorporate this condition to design Q 0 (z) [7]. The condition is derived below. 

Let N r , N p , and N q be, respectively, the orders of Rq(z), Pq(z), and Qo{z ). Let the impulse responses 
of Rq(z) ) Po(z), and Qo(z) be ro (n), po(n ) and <3 , o( n ) ) respectively. Then r*o(n) is the convolution of po(^) 
and q 0 (n), i.e., 


r 0 {n) - n — m)q 0 (m) (21) 

m=0 


Given p 0 (n ), we want to constrain q 0 (n ) so that r 0 (n) is a Nyquist (2 M) filter. In this case, q 0 (n) must 
satisfy 


Y. Po(2Mn - m)q 0 (m) = c6(n - n 0 ) (22) 

m=0 


for some n 0 and c. Constant c is a scalar factor; we can consider c = 1. The equivalent of Eq. (22) in 
matrix form is 


r = Pq (23) 

where r is a ([(A^, -f- N q ) /2M J + 1) x 1 column vector with the noth entry equal to 1 and all others equal 
to 0, P is a ([ (N p + N q )/2M\ + 1) x (N q + 1) matrix with 

rpi _ fpo(2Mm - n + 1), if2Mm-n+l>0 
mn \ 0, otherwise 


and q = [q 0 ( 0) ? 0 (1) • • • g 0 (iV g )] T - 

The condition in Eq. (23) can be easily incorporated in the eigenfilter approach [7]. 

B. A Shortcut Design 

Let Pq{z) be a Nyquist (2 M) lowpass filter with zero phase. Let Qq(z) also be zero phase with its 
passband covering Pq(z ) as depicted in Fig. 15, i.e., Qo(e^ u> ) ~ 1 in the passband and the transition band 
of P 0 (e JW ). Because R 0 (z) differs from P 0 (z) only in the stop band where the magnitude of P 0 (z) is small, 
we have 
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Fig. 15. The passband of Qq(z) covers the passband and transition band of Pg(z). 


R 0 (z) « P 0 (z) 


( 24 ) 


In this case, Ro(z ) will still be close to a Nyquist (2M) lowpass filter. 

In the above discussion, Pq{z ) and Qq(z) are zero phase and, hence, noncausal. Some delays can be 
added to make the filter bank causal, since both prototypes are FIR. In the two following methods, the 
synthesis prototype Qq(z) is chosen as the time-reversed version of P 0 (z), i.e., Q 0 (z) = z~ N r P n (z) where 
N p is the order of Pq{z). 

C. Nonlinear Optimization of P Q (z) 

The filter coefficients of P 0 (z) can be directly optimized to minimize the stop-band energy of P 0 (z) 
subject to the Nyquist condition in Eq. (22). The stop-band energy of P 0 (z) is 


<ps= r 

J(n/2M)+e 


\Po (e JW ) \ 2 dw 


But with Qo(^) — ^ Np Po(z), the right-hand side of Eq. (22) becomes a quadratic form of po(n). To 
minimize <f> a , we need a nonlinear optimization package that can incorporate nonlinear constraints. The 
detail of this design technique is documented in [12]. 

Instead of optimizing the coefficients of Pq(z) with nonlinear constraints, we can use a objective 
function to reflect how close Ro(^) is to a Nyquist (2 M) filter. A possible objective function is 

pn/M _ 

<f>p = J q (|Po (e j “) | 2 + |P 0 (ef(«-(*/"))) j 2 - i) dw 
Combining both (p s and <j> p , let 


0 = a<j)p + (1 - a)cj) s 

where 0 < a < 1. The new objective function 0 can be minimized using the usual nonlinear optimization 
programs (e.g., [13]) without any constraints. 

D. Kaiser Window Design 

In this approach, P 0 (z) is a lowpass filter with an even order obtained through Kaiser window design. 
That is, 


Po{n) 


sin(o ) c n) 
nn 


w(n ) 
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where w(ri) is a Kaiser window [7]. After we choose the stop-band attenuation and the width of the 
transition band, the length of the window can be estimated by a formula developed by Kaiser. In this 
case, the window is completely determined [7]. The cutoff frequency u) c of the ideal filter is the only free 
parameter left for the design of Pq(z). 

According to Eq. (22), we choose a simple objective function, (p Kaiser) to reflect the closeness of Rniz ) 
to a Nyquist (2 M) filter: 


< pKaiser = max \r 0 (2Mn ) - din - n 0 )| 

Experiments show that (j) Kaiser is a convex function of w c ; we can adjust the parameter cu c to find the 
best P 0 (z), which yields the smallest </> Kaiser- Design examples will be given in the next section. 

E. Comparison of the Four Design Methods 

Of all the four approaches stated above, the shortcut design and the Kaiser window design are probably 
the easiest. In the shortcut design approach, we only need to design two lowpass filters with certain 
specifications described above. In the Kaiser window design, only the parameter w c needs to be optimized, 
yet no nonlinear optimization is involved. Also, the value of the objective function (j) K aiser at a given w’ 
can be computed easily. 


V. Design Example 

We now present one design example of the under-decimated system. The type 1 cosine modulated filter 
bank is used in this example. The Kaiser window approach is adopted for the design of the prototype 
filter. 


Example 1: Tunable Multilevel Filter. The system has 20 channels. In this case, M — 10. The 
analysis bank prototype filter P 0 (z) is linear phase with order N p = 120, stop-band attenuation 85 dB, 
passband edge co p — 0.047T, and stop-band edge uj s — 0.0997T. The synthesis bank prototype Qo(z) is 
chosen as the time-reversed version of P 0 (z). As elaborated in Section III, the resulting overall response 
will have linear phase. Figure 16(a) shows the magnitude response of Pq(z). The normalized magnitude 
responses of the analysis filters are shown in Fig. 16(b) and (c). The synthesis filters are time-reversed 
versions of the analysis filters; the magnitude responses of the synthesis filters are the same as those of 
analysis filters. 

After designing the prototype filters, the values of a*, can be changed freely to obtain the desired 
overall response, T(z). For instance, we set a 0 = = 1, a 2 = a 3 = a 4 = 0, a 5 = a 6 = a 7 = 0.7, and 

a s =a 9 = aio = 0.3. The magnitude response of the resulting T(z) is plotted in Fig. 16(b). Since T(z) 
has linear phase, we did not show the phase response. The corresponding dB plot of Fig. 16(d) is shown 
in Fig. 16(e). 
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Fig. 16. Example: (a) the magnitude response of the prototype filter, Pq(z), (b) the normalized magnitude 
responses of the analysis filters in the first subsystem, (c) the normalized magnitude responses of the 
analysis filters in the second subsystem, (d) the magnitude response of the overall response T{z), and (e) 
the magnitude response of the overall response Rz) in a dB plot. 
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Appendix 

Implementation of the Type 1 Cosine Modulated Filter Bank 


In this appendix, we prove that the implementation cost of the analysis bank of the type 1 cosine 
modulated filter bank is that of the analysis prototype filter plus two DCT matrices. 


Let 


Po(z) 


2M — 1 


E G - C M ) *- 

n=0 


where G n {z) is the nth type 1 polyphase component of P 0 (z). Then 
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(A-l) 


2M-1 

Pk(z) = Gn ( z2M ) Z~ n W~ kn 

n = 0 

Rewriting analysis filters in Eq. (13) in terms of the polyphase components of Pq{z) with a k ,b k , 
b' k as in Eq. (15), we obtain 


2M—1 

H k (z) = 2 a k G n (z 2M ) 


n—0 


2M-1 

H’ k {z) = 2 £ a k G n (z 2M ) 


71 = 0 



k = 0,1, ■ • • , M 


k = 1, 2, • • ■ , M - 1 


Define a 2M-component vector h (z) given by 


( H q (z) \ 


h (*) = 


H m {z) 

H[(z) 


Wm-i(z)J 


Using Eq. (A-2), h(z) can be written as 


h (z) = 2 



OUC A 0 C\ fg 0 (z 2M ) 
Ai ) \ S A X S ) \ 0 


0 \( e{z) \ 

gi (z 2M )j[z- M e(z)J 


where g i(z), A i, and A* are diagonal matrices with 


[go(*)]fcfc = G k {z), [gi(.z)3fcfc = G k+M (z), k = 0, 1, • ■ • , M — 1 


[Aojfcfc = (-l) fc , k — 0,1,- • ■ ,M 
[Ai]fcfc = (-l) fc , fc = 1,2, • • • ,M — 1 

[Aojfefc — a k , k = Q,l,---,M 

[Aijfcfc = c i k , k = 1, 2, • • • , M — 1 
and C and S are (M + 1) x M and ( M — 1) x M matrices with 


a ' k , and 


(A-2) 


(A-3) 


(A-4) 


(A- 5) 


(A-6) 
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7T 

[C] mn = cos(— mn), m = 0, • ■ • ,M, n = 0, • • • ,M - 1 


[S]mn = sin(— mn), 


ra = 1, • • • , M — 1, n = 0, • • • , M — 1 


(A-7) 


Define two M x 2M matrices 


T 0 = (C A 0 C) and Ti = (S AiS) 

Prom Eq. (A-3), we can draw Fig. A-l, a schematic implementation of the 2 M channel cosine modulated 
system. The input to To and Ti, a(n), can be partitioned into two M x 1 vectors: 


a(n) = 


( a o(n)\ 

U (n)J 


Their dependence on n will be dropped for convenience. As indicated in Fig. A-l, d 0 and di are the 
outputs of T 0 and T 1 , respectively: 


d 0 = T 0 a and dj. = Tia 
From the definitions of T 0 and Ti, we know 

do — Cao + AoCai and di = Sao + AiSax 
In [14], it is pointed out that C and S have the following properties: 


(A-8) 



Fig. A-1. A schematic for the implementation of the analysis bank of the 
type 1 under-decimated cosine modulated 2Af-channel filter bank. 
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A ° c = c (o jL) + ^ 0 • ») - d AlS = - S (o jZJ 

Using these relations, Eq. (A-8) becomes 

d „ = C (a„+(o JM. 1 ) ai ) + [ai,or ' d i = s(^>-(o (A-9) 

where [ai] 0 is the first element of ai, and r = [1 - 1 ■ ■ • jt Equation (A-9) allows us to have a 

more clear idea of the implementation of T 0 and Ti . The implementation of the synthesis bank is similar 
to the above. 

Now we can draw a more detailed and complete picture for the implementation of the 2 M channel 
filter bank. For simplicity, we can choose the synthesis filters to be the time-reversed versions of the 
corresponding analysis filters without the scalars < 2 The efficient implementation of the filter bank is 
shown in Fig. A-2. From Fig. A-2, we can see that the major complexity of the analysis bank is the cost 
of the analysis prototype plus matrix C and matrix S. Matrices C and S can be implemented by fast 
algorithms for DCT and DST matrices [11]. The process is similar for the synthesis bank. 



C (M+1) x M 


S (M-l)xM 

a 0 , a 1 square matrices 

Fig. A-2. Efficient implementation of the type 1 under-decimated cosine modulated 2/W-channel filter bank, 

where k= N- 2M + 1 and e 0 = [10... 0]^. 
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